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Abstract 

Various fluid mechanical systems enjoy a hidden, higher-dimensional dynam- 
ical Poincare symmetry, which arises owing to their descent from a Nambu- 
Goto action. Also, for the same reason, there are equivalence transformations 
between different models. These interconnections, summarized by the diagram 
below, are discussed in our paper. 
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I. INTRODUCTION 

In this paper we shall be concerned with nonlinear dynamical systems that are described by 
a density of "matter" p, flowing in time {£} with local velocity v on a rf-dimensional surface 
coordinated by {r}. The vectorial nature of v is not unrestricted: v is a function of V0, 
where 9 is a velocity "potential" , and we shall examine several such functions. (When v is 
linear in V0, the flow is irrotational, V x v = 0.) The density is linked to the velocity by a 
continuity equation involving the matter current j = vp, 

p(£,r) + V-(v(£,r)p(£,r))=0 (1.1) 

while the velocity satisfies an "Euler" equation 

v(£, r) + v(£, r) • Vv(£, r) = f (p, v) . (1.2) 

[The over-dot always indicates differentiation with respect to the temporal argument, while 
the gradient V (unless further specified) differentiates the spatial argument.] 

We shall examine theories with various expressions for the force (per unit mass) f , which 
lead to Galileo, Poincare and additional unexpected kinematical symmetries of the equations, 
and which sometimes produce completely integrable systems, with an infinite number of 
local conserved quantities. The existence of a velocity potential 6 allows the above equations 
to be formulated with an action principle, which is usually unavailable for this purpose. 
Consequently the symmetries that we find are in fact Noether symmetries, which leave the 
action invariant. 

Additionally, we shall present limiting and equivalence transformations between different 
models, which allow mapping solutions of one model onto solutions of another. 

Various topics that we discuss have already appeared in the literature. A common fea- 
ture unites the diverse models that we study: they have a common antecedent in that 
they can be gotten from a parametrization-invariant Nambu-Goto action for a d-brane on 
a d + 1-dimensional space, moving in (d + 1, l)-dimensional space-time. [A "d-brane" is a 
(i-dimensional extended object: d = 1 is a string, d = 2 is a membrane, etc. A d-brane inhab- 
iting d + 1-dimensional space divides that space in two.] When a light-cone parametrization 
is selected for the Nambu-Goto problem, one derives the Euler and continuity equations 
for a d-dimensional "Chaplygin gas" [in Eq. (|1 . 2| ) , f oc -^Vp]. Alternatively, a Cartesian 
parametrization produces the (i-dimensional "Born-Infeld" model (see below). 

The relation between membranes and planar fluid mechanics (the d = 2 case) was known 
to Goldstone 0, and was developed by his student Hoppe (with collaboration) 0. Similar 
connections, yielding equations in one spatial dimension, were discussed by Nutku (and 
collaborators) ||. Here we consider the general (i-dimensional case, and use the common 
ancestry of the various fluid-mechanical models to posit unexpected transformations between 
them, and to identify hidden, dynamical symmetries in each model, which derive from the 
high degree of symmetry of Nambu-Goto parent theory. 

For d = 1, the models that we study become especially simple for two reasons. First, their 
common antecedent is a string (1-brane) moving on a plane (2-space), and for this system 
the Nambu-Goto equations are integrable J3J. Second, the requirement that velocity be 
expressible in terms of a potential poses no restriction in one dimension, where any function 
can be related to the derivative of another function. In this way one may understand that 
the d = 1 models are completely integrable, as has been noted previously 



In Section II we consider noninteracting systems, with Galileo- and Poincare-invariant 
kinetic terms. Specific interactions that preserve Galileo and Poincare symmetry, as well 
as higher, dynamical symmetries are discussed in Section III. The relation of these to the 
Nambu-Goto theory is explained in Section IV, where we also exhibit mappings between 
the Galileo-invariant and the Poincare-invariant models. The last Section V is devoted to 
models in one spatial dimension. 



II. FORCE-FREE MOTION 



The force-free problem, f = 0, describes the free flow of dust. Eqs. ( |1 . 1| ) and ( |1.2| ) are readily 
solved in terms of initial data, specified (without loss of generality) at t — 



v (t, r ; 



i=0 



i=0 



Po(r) 



v oir J 



Upon defining the retarded position q(t, r) by the equation 

q + tv (q) = r 
one verifies that ( |l.lj ) and ( |1.2| ), with vanishing right side, are solved by 

v(t,r) = v (q) 



p(t,r) =p (q) 



det 



dq l 



(2.1) 
(2.2) 

(2.3) 

(2.4) 
(2.5) 



The free equations, with v restricted to a function of V#, possess a variational action 
formulation, which was first given by Eckart for Galileo-invariant nonrelativistic motion |j]. 
We reproduce and generalize his argument. 

The Lagrangian for N point-particles of mass m in free nonrelativistic motion is the 
Galileo- invariant kinetic energy ^J2n=i mv n{t)- ^ n a continuum description, the particle 
counting index n becomes the continuous variable r, and the particles are distributed with 
density p, so that J2n=i v n(t) becomes J d d r p(t, r)w 2 (t, r). But we also wish to link the 
density with the current by the continuity equation (|1 . 1|) , which can be enforced with the 
help of a Lagrange multiplier 9. We thus arrive at the free continuum Lagrangian 



L 



Galileo 



d d r 



p\mv 2 + 6(p+ V- (vp) 



(2.6) 



Since L is first-order in time, and the canonical 1-form J d r 9pdt does not contain v, v may 
be varied, evaluated and eliminated [[/[. We find 

pmv - pV6 = (2.7) 

showing that V# is the local momentum p = mv, and the velocity is irrotational: 

v= — V6 

m 



V x v = . 
3 



{21 



The Lagrangian (2.6) becomes 



j Galileo 
L 



dV 



- p 



2m 



(2.9) 



where the subscript denotes absence of interaction. 

Varying 9 in ( |2.9j ) regains the continuity equation ( |1 . 1| ) , while varying p produces the free 
"Bernoulli" equation for the velocity potential 9 



«+2^ V 



. 



(2.10) 



This is also recognized as the free Hamilton- Jacobi equation. The gradient of (|2.10| ) gives 



rise, in view of (2^), to the free Euler equation ( |1.2|) (with f = 0). 

Remarkably the same equations emerge for a kinetic energy T that is an arbitrary function 
of v. If we generalize ( |2.6| ) to 



L Q = J d<V [pT(v) + 9(p + V- (vp) 



we get, instead of (|2] 



<9T(v) 
<9v 



p = V9 



and (I2.9T) becomes 



d d r 



/Mv-^^-T(v 



(2.11) 



(2.12) 



(2.13) 



where it is understood that the Legrendre transform of T is expressed in terms of V# by 
inverting fl2.12|) . 



Varying 9 in ( |2.13j ) again gives the continuity equation, 

Mo . f , d __ 6 (dTiy) 







p — d d r pv 



59 r J " ' r ' 69 \ dv 

= p + V • (vp) 
while varying p leaves a generalization of the free Bernoulli equation: 

SL 







5p 



dT 
-0-v.-+T ( v) 



(2.14) 



(2.15) 



With the help of ( p.!2| ), this implies 



P 

or 1 



, d 2 T 

, d 2 T 
drWv 1 



, d 2 9 

, dv k d 2 T 



g r j dv k dv l 



(2.16a) 



On the other hand, from ([2.1 2|) it follows that 



d ■ d 2 T 
dr l dv l dv k 



Eqs. ( |2.16a|) and (|2.16b| ) are consistent only if the free Euler equation (|1.2|) holds (provided 



the matrix » a T ■ has an inverse 

av l av-> 



With a general form for T(v), the local momentum (2.12) p = |- remains irrotational 



while the velocity, as determined by inverting ( |2.12|) , becomes a nonlinear function of V0. 

Evidently the solution ( |2.1| )- fl2~5| ) works with arbitrary kinetic energy, whose specific 
form enters only in the fixing the relation between v and V(9. However, the initial data for 
the velocity must be consistent with the expression of the velocity in terms of V0. 

One may present a family of constants of motion: 



C 



|d d rp(t,r)C(v(t,r),r-tv(t,r)) . (2.17) 



The time independence of C is established either by differentiating with respect to t and 
using the free equations of motion to prove that ^ = 0, or more easily, by inserting the 
solution ( |2.3| )-( |2~3| ) into ( |2.17| ) and changing integration variables from r to q. (Carrying out 
these manipulations requires assuming that po and vo obey appropriate regularity conditions 
and drop off sufficiently at large distances.) 

Various constants of motion arise from invariance against time and space translation (en- 
ergy E and total momentum P, respectively) as well as space rotation (angular momentum 
Lij), provided T(v) carries no explicit time and coordinate dependence, and does not depend 
on any external vectors, i.e., T(v) = T(v). These constants are 

E = H= I d d rU H = p(v ■ — - T(v)) (2.18) 

P = / d d r V V = pV6 = p^-t = pp (2.19) 

J <7V 

L id = f d d r {r ir P j - r^V 1 ) . (2.20) 

Also shifting 9 by a constant is a symmetry, leading to conservation of 

N=fd d rp. (2.21) 

To recognize that these constants of motion involve particular forms for C(v, r — tv) in 
( p.!7| ), we recall that according to ([2.12|) V6 1 is a function of v, and the two are colinear 



when T(v) = T(y). 

The densities Ti and "P are components of an energy-momentum tensor, T 00 and T° l 
respectively, which satisfy continuity equations with energy flux T*° and momentum flux T* J 

T 00 = H (2.22a) 

T i0 = v l H (2.22b) 

T° l = V 1 (2.22c) 

T ij = v i<pj _ (2.22d) 

[T 2J is symmetric when T(v) = T(v).] The continuity equations 
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d 
T 00 + — T l ° = (2.23a) 

dr l y ' 

d 
T° 3 + — -T 3 = (2.23b) 

or 1 

are entirely equivalent to the free dynamical equations ( |1 . 1| ) and ( |1.2|) , with vanishing force. 
The symplectic structure, which is determined by the canonical 1-form in (|2.13|) , indicates 
that the only nonvanishing bracket is J7| 

{e{t,v) )P {t,v')} = 8{v-v') (2.24) 

or equivalently 

{p(t,r),p(t,r')} = V8(r-r>). (2.25) 



With these, one verifies that the constants of motion (|2.18| )-( [2.21| ) generate the appropriate 
infinitesimal transformation on the variables 9 and p. 

Specific forms for T(v ) support additional, kinematical symmetries and lead to further 
constants of motion. In the nonrelativistic case presented in Eqs. ( p.6| )-( pTT0| ), we have 
Galileo invariance against boosts by velocity u. The transformation law for the fields reads 

p(t,r) -> p u (t,r) = p(T,R) 

6{t, r) -► 9 u (t, r) = 6{T, R) + m(u ■ r - u 2 t/2) (2.26) 

where the transformed coordinates are boosted: 

t^T = t 
r ->R(i,r) = r - tu . (2.27) 

The inhomogenous terms in 9 U are recognized as the well-known 1-cocycle of field theoretic 
realizations of the Galileo group. Also they ensure that the transformation for the velocity 
v = V# jra 

v(t, r) -> v u (t, r) = v(t, r - tu) + u (2.28) 

is appropriate for the co-moving velocity. 

The conserved quantity arising from the Galileo symmetry is 



B 



iP — m I d d r rp 

-m fd d r(r-tv)p (2.29) 



where the last equality casts B in the form ( ^.17[ ). With the help of the bracket (2.24), B 



generates the infinitesimal transformation on the fields, and its bracket with P closes on N, 
thereby exposing the familiar Galileo 2-cocycle, which provides an extension of the algebra: 

{B\P j } = 5 ij mN . (2.30) 

The free Galileo-invariant theory possesses further symmetries, which survive even in the 
presence of a particular interaction. Hence we postpone discussing them until later, when 
interactions are included. 
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In the subsequent, in addition to the Galileo-invariant case, we shall also be concerned 
with a relativistic, Poincare-invariant model for which the point-particle kinetic energy is 
—mc 2 J2n=i \A — v n(t)/ c2 - Upon passing to a continuum description, as in the nonrelativistic 
case, we find 

T(v) = -mc 2 ^l-v 2 /c 2 (2.31) 

dT(v) mv „ . 

— ^ = p = = V6 2.32 

<9v ^1 - v 2 /c 2 



eve 



m 2 c 2 + (vey 



leading to 



mc 2 



and Lagrangian 



L Lorcntz = j ^ 



^(T,R) + Cv /mV + (V^(T,R)) 
p(t,r) -* p u (t,r) = p(T,R 



:» 



\2 



with Lorentz-boosted coordinates 



t->T(t,r) = £cosh/3 + -/3-rsinh/? 

c 



(2.33) 



W = pcJm 2 c 2 + (V0) 2 = pcJm 2 c 2 + p 2 = p , (2.34) 

^1 - w 2 /c 2 



(2.35) 



- pc^m 2 c 2 + (V6») 
In the nonrelativistic limit this becomes 

^Lorcntz _> _ mc 2,y + ^Galileo _ ^gg) 

Under Lorentz boosts, with velocity u, the fields transform as 



d t 9(T, R) + c^m 2 c 2 + (V r 0(T, R))< 
6(t,r)^6 u (t,r)=6(T,R) (2.37) 



r ->R(t,r) =r + 3(ctsinh/5 + 3-r(cosh/3-l)) (2.38) 

where (3 = u/c. Invariance is most easily verified by writing the action corresponding to 
( p| ) as 

jLorentz = _ f dt ^ ^ + C J^2 + ( V fl)2 | . ( 2 .39) 



The infinitesimal version of the field transformation (|2.37| ) 



^ = /3 .(__ + rf V)p-,3.-p 

se = P'(- — + ctv)e (2.40) 
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is generated by the Lorentz constant of motion 

L = tP - f d d r rH/c 2 



f d d r (tpVO - -p^Jm 2 c 2 + ( V0) 

-m / d d r (r - tv)p = (2.41) 

^ J\ — v 2 /c 2 



with the last equality exhibiting the form ( |2.17|) . The transformation laws for p and 9 



ensure that j a = (l,v/c)p and U a = (1, v/c) , = = transform as Lorentz vectors, so 



that pJl — v 2 /c 2 and 6 are scalars ||. The equation of motion ( j2.15|) for together with 
( |2.31| )- (|2.33| ) implies that the Lorentz vector d a 6 satisfies a Lorentz-invariant equation 

{djf = m 2 c 2 . (2.42) 



Note that the Lorentz transformation law ( |2.37|) for 9 does not involve a 1-cocycle, which 



is a nonrelativistic effect. It is interesting to see how (|2.26|) arises in the nonrelativistic 



limit. By comparing the relativistic action ( |2.39| ) to the nonrelativistic one, we see that the 
relationship between #r and #nr, — the relativistic and nonrelativistic variables — is 

#r = #nr - mcH . (2.43) 

Applying the Lorentz transformation law ( |2.37| ) to 6^ implies that ^nr(^j r ) — mc 2 t —> 
6*nr(7\ R) — mc 2 T or 6*nr(£, r) — > Onr(T, R) + mc 2 {t — T). We evaluate the nonrelativistic 
limit of the last quantity from (|2.38|) and find mc 2 (t — T) — > m(u • r — u 2 t/2), which matches 
the 1-cocycle in (|2.26f) . 



Similar to the Galileo-invariant theory, this Poincare-invariant model possesses further 
symmetries, which we shall discuss below, when we include an interaction that preserves 
them. 



III. MOTION IN THE PRESENCE OF INTERACTIONS 

A. Nonrelativistic motion 

Interactions that preserve the Galileo symmetry of the free, nonrelativistic motion can be 
included by adding a ^-independent potential V(p) to the Lagrangian (|2.9| ): 



L Galileo = I d d r ■• «'V^- 



J V 



P^t-v(p) 



(3.i; 



With nonvanishing V, ( p.3| )-( p3|) are no longer solutions, and the quantity ( 2.1 7|) with 



arbitrary C(v, r — tv) is no longer constant. Of course the Galileo generators ( |2.18|) - (|2.21 
and ( PJSD , with 



remain time independent. The energy-momentum tensor retains the form (|2.22|) , with 7i 
given by (|3.2|); in T i0 , H of ( |3.2| ) is diminished by V — pfr^; T° l is unchanged and T u 

acquires the addition 5^7 V — pfp)- 

The dynamics implied by (|3.1|) arise in diverse physical contexts. The most directly 
physical application is to isentropic, irrotational motion in fluid mechanics with the "force" 

f(p) = t^ corresponding to the enthalpy and (p a 2 ) is the speed of sound JJ. Al- 
ternatively one finds ( |3.1| ) (with V depending also on Vp) in the hydrodynamical formulation 
of quantum mechanics, which emerges when the wave function is presented as (l(| 



^ = pV2 e ^ (33) 



[In this context, the inhomogenous Galileo transformation of 9 ( 2.2(f ) corresponds to the fa- 



miliar change of phase in a wave function under Galileo boosts, while shifting 6 by a constant 
is just the phase-invariance of quantum mechanics, which leads to probability conservation, 
i.e., constant N in ( |2.21| ).1 

But we shall be especially concerned with the case 

V(p) = \/p (3.4) 

which arises in the study of "d-branes" - d- dimensional extended objects - - moving on 
a (d + l)-dimensional space, in (d + 1, l)-dimensional space-time, and descending from a 
Nambu-Goto action (see Section IV) |I]j@]. Furthermore, Q3.4J ) arises in the nonrelativistic 
limit of a Poincare-invariant model with interactions, which we shall also describe below. 
The equations of this theory, which follow from 

^Galileo = f d d r [ Q - p _ (yOf _ A 1 (35) 

J I 2m p) 



f> -V.(— p)---() (3.6) 



read in their Bernoulli form 



rn 

« + ^^4 P.7) 

2m p" 1 

while their Euler form is gotten by recalling that v = V^/m 

p + V • (vp) = (3.8) 

— 2A 

v + v-Vv = -Vp . (3.9) 

mp^ 

These are the equations for a "Chaplygin gas." 

For this model there exist further symmetry transformations [11]. The action jG allleo — 
/ d£L^ allleo is invariant against a rescaling of time, with parameter u, under which the fields 
change according to 

p(t,r)^p u (t,r) = e- u p(e u t,r) 

6{t, r) -> u (t, r) = e^O^t, r) (3.10) 

9 



and the time-independent generator of the infinitesimal transformation is 

D = tH - f d d rp9 . 



(3.11) 



Furthermore, the action is also invariant against a field-dependent diffeomorphism, implicitly 
defined by 



t -> T(t, r) = t + u ■ r + ±oj 2 9(T, R)/m 
r -► R(£, r) = r + u>0(T, R)/m 

where the transformed fields are 



(3.12) 



p(t,r) -* p u (t,r) = p(T,R)- 
0(*,r)->0 w (t,r) = 0(T,R) 

and |J| is the Jacobian of the transformation. 



J I 



(3.13) 



J = det 






1 - — ■ V0(T, R) - — 0(T, R) x 
m 2m 



(3.14) 



Here a; is a (vectorial) parameter of the transformation, with dimension of (velocity) l , and 
the time-independent quantity 



G = f d d r (vH - eV/m) 



(3.15) 



generates the infinitesimal transformation. 

Note that the generators D and G remain time-independent even in absence of the 
interaction ( p.4| ), hence these symmetries are also present for the free theory. The generators 
are not of the form ( |2.17| ): they involve 6, and cannot be written in terms of v = V0/m. 
Finally we note that bracketing of the additional generators with the Galileo generators 
of the nonrelativistic theory on a (d, l)-dimensional space-time produces an algebra which 
is isomorphic to the Poincare group in (d + 1, 1) dimensions, @,|l^] under which (t, 9, r) 
transforms as a (d + 2)-Lorentz vector X^ in light-cone components (t = X + , 9 = X~) M. 

Using (|3.7| ), we may eliminate p, and describe the model solely in terms of 9, whose 
dynamics is governed by the Lagrangian 



-2VX 



d d r\ 9+ (V 



2m 



(3.16) 



It is seen that the "interaction strength" A in fact disappears from the equations of motion for 
9; A serves only to normalize the Lagrangian. In the free theory it is not possible to achieve 
this compact formulation. Furthermore, the dynamical equations can be summarized by an 
equation for 9, which follows from ( |3.6| ), once ( |3.7| ) is used to eliminate p, or alternatively 



the equation is derived from (ft.l6|) 



Of 



( 



V- 



9 + 



V9/m 



\ 



2m 
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(3.17) 



In spite of their awkward appearance, ( p,16|) and ( |3.17|) are Galileo invariant in (d, 1) space 



time, and possess a hidden, nonlinearly realized Poincare symmetry in (d+ 1, 1) space-time 
(which is a descendant of the symmetries of the Nambu-Goto action; see Section IV). 

Apart from the intrinsic interest in this nonlinear realization of a kinematical/dynamical 
Poincare symmetry, which is provided by (|3.10|) - (|3.15|) supplementing the linearly realized 
Galileo symmetry, the new symmetries have the useful consequence of generating new solu- 
tions to the equations of motion ( |3.6| )-( |3~9| ) from old ones. For example, the time-rescaling 
invariant, particular solution 

0{t,r) = 



mr 2 



2(d-l)t 

2A 

P(t^) = \l—Ad- 1 )- (3-18) 




can be transformed by ( |3.12| )-( |3.14| ) into new solutions 9 U and p u , which are very different 



in character from (|3.18|) ||1 1|| . Note that in ( p,18|) we must have d > 1 and A > 0. 



At d = 1, we can obtain general, time-rescaling invariant solutions. With the Ansatz 
9 oc 1/t, ( |3.17|) leads to a second order differential equation for the x-dependence of 9. 



Therefore solutions involve two arbitrary constants, one of which fixes the origin of x, and 
can be ignored. The other, which we call k, appears in two distinct families of solutions 
(which are related by an imaginary shift of x) : 



m , 2 



9(t,x) = — cosh kx (3.19a) 

2k 2 t 



m . , 2 



6(t, x) = —— sinh 2 kx . (3.19b) 

ZsKj l 



For real 9, k must be real or imaginary. When a real p is computed from ( |3.7|) , we find that 
k must be real for A > 0, imaginary for A < 0, and a nonsingular density exists only with 
( Ey]§ ) for A > 



^^VifJjte- (3 - 20) 

The current j = —p exhibits a kink profile (derivation with respect to the single spatial 
variable is indicated by a dash) 

/2A 
j(t,x) = T\ — tanhkc (3.21) 

V m 

where the sign is fixed by the sign of t. 

In the last section we shall further review the d = 1 case. 

Another interesting solution, which is essentially one- dimensional, even though it exists 
in arbitrary spatial dimension, is given by 

0(t, r) = 6(fi • r) + mil • r - \mt{u 2 - (n • u) 2 ) . (3.22) 

Here n is a spatial unit vector, and u is an arbitrary vector with dimension of velocity, 
while G is an arbitrary function with static argument, which can be boosted by the Galileo 
transformation (|2.26|) . The corresponding charge density is time-independent: 
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j2X/m 

P{t f r) = - 1 (3.23) 

n ■ u + B (n • r)/m 



and the static current becomes 



2X / u — fi(fi • u) 

j(t,r =V — ( fi +5 , n/ V ^ • (3-24) 

V m V n • u + B '(n • r)/m; 



B. Relativistic motion 

We seek an interacting generalization of LQ° rentz , which preserves Poincare invariance. To 
find this, proceed as follows. Let 



Lorentz 



Jd d r{9p-H(p,p)} (3.25) 



with p given by fl2.32| ) and 7i to be determined. The symplectic structure is as in the free 



theory, hence the Poisson brackets retain the form ( |2.24|) , ( |2.25| ). We calculate the Poisson 
bracket between two Hamiltonian densities; in one the fields are evaluated at r, in the other 
at r': 

{H(r),H(r')} = |dr''dr™{^ • W(r" - r"')g^j - r ~ r] . (3.26a) 

We assume that 7i is a local function of p and p, so that the functional derivatives lead to 
ordinary derivatives, and ( |3.26a|) becomes 

On the other hand, the Dirac-Schwinger condition for Lorentz invariance states that the 
bracket ( |3.26|) should give rise to c 2 times the momentum density T 5 , which in this problem 
is given in Q2.19| ) as 

V = pp . (3.27) 

Rotational invariance requires that the p dependence of TC is only on the magnitude p. Thus 
we conclude that 

4 fiN^ c2 - (3 - 28) 

While many forms for Ti can solve ( P-28| ), we take a solution that is relevant to the present 
context, i.e., it generalizes in a simple manner the free Hamiltonian density Q2.34J ), leads to a 
theory that descends from the Nambu-Goto action, and coincides in the nonrelativistic limit 
with the Galileo-invariant Chaplygin gas model: 

2 

H = ^Jp^c 2 + a 2 ^m 2 c 2 + (V9) 2 = y/ ' p 1 + a 2 /c 2 m ° = . (3.29) 

J\ — v 2 /c 2 

12 



Evidently the parameter a measures strength of "interaction". 

Thus an interacting, Poincare-invariant theory is described by the Lagrangian 



■ Lorentz 



d d r 



p 2 c 2 + a 2 Jm 2 c 2 + (V6 



(3.30) 



The corresponding conserved Lorentz generator takes the same form as in the first equality 
of fl2.4ip , with 7i given by ( |3.29| ), and it generates the infinitesimal transformation 



r d 



Sp = (3-(-- + ctv)p - (5 ■ -^p 2 + a 2 /c 2 



■cdt 



(3.31) 



The finite transformation law is gotten by iterating (|3.31|) , and the generalization of ( p. 37 ) 
becomes 



pit, r) - p u {t, r) = p(T, R)|(fi + + «-) + \JP 2 (T, R) + a 2 /c 2 ±(ft+ - fi_) 
6(t,r)^6 u (t,r) = 9(T,R) (3.32) 



where the Lorentz transformed coordinates (T, R) are given in ( 2.38 ) and 



0(T, R) ± c^Jm 2 c 2 + (V0(T, R)) 2 
d t 6(T, R) ± Cv /m 2 c 2 + (V r 0(T,R)) 2 

It follows that j a = (p, ^p 2 + a 2 /c 2 ) and U a = ( 

while 9 is a scalar. 

The equations of motion are 

eve 



-yV+aVc 2 ' c / a/1-^/c 2 



(3.33) 



are Lorentz vectors, 



p + V 



m 2 c 2 + (V#) : 



^p 2 + a 2 /c 2 ) = p + V • (vyp 2 "^ 2 /^ 2 ) = 



'm 2 c 2 + (V#) 2 . 
+ pc " ; =- = + 

p 2 + a 2 /c 2 



?71C 



p 2 + a 2 /c 2 w 1 — t> 2 /c 2 



(3.34a) 
(3.34b) 



Using ( 3.34b ) to express p in terms of 6, 

a d Q 9 



P 



c 2 ^m 2 c 2 - (d^Of 



(3.35) 



and substituting this in ( |3.34aj) yields a second order, Lorentz covariant equation for 9. That 
equation may also be gotten by eliminating p from ( ^.3U| ) and deriving a Lagrangian for 9. 



L e = -a f d d r sjm?c 2 - {d a 9f 



This is the "Born-Infeld" Lagrangian, leading to the equation of motion 

1 



d c 



^m 2 c 2 - (d^9) 2 
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d n 9 ) = . 



(3.36) 



(3.37) 



As in the nonrelativistic theory [see (|3.16|) ] the possibility of expressing p in terms of 



9 requires presence of the "interaction," a ^ 0, whose nonvanishing strength disappears 
from the nonlinear, interacting equations for 9. The energy-momentum tensor for the theory 
( |3.3U| ) is Lorentz covariant, of second rank, and symmetric. After eliminating p with ( |3.35| ), 
the resulting expression, depending solely on 9, bears the usual relation to its Lagrangian 

(H). 

Since the interacting Lorentz-invariant model is a descendant of the Nambu-Goto La- 
grangian (see Section IV), it comes as no surprise that it too possesses additional kinematic 
symmetries, whose generators supplement the generators of the linearly realized Poincare 
group in (d, 1) dimensions to give a nonlinear realization of dynamical Poincare algebra in 
(d + 1,1) dimensions [|H| . 



The additional symmetry transformations, which leave (|3.30|) or ( |3.36| ) invariant, involve 
a field-dependent reparametrization of time, defined implicitly by 

*->T(t,r) = — - t - + ^ y tanh mc 2 uj (3.38) 

cosh mc 2 u mc z 

under which the field transforms according to 

9(T r) 

6(t,r) -► 6Jt,r) = — y -^ mc 2 t tanh mc 2 uo . (3.39) 

cosh mc l ijj 

[We record only the action of the transformations on 9; their effect on p can be read off from 
( p.35| ).1 The infinitesimal generator, which is time independent by virtue of the equation of 



motion (|3.37|), is 



D = f d d r (m 2 cHp + 9^/p 2 c 2 + a 2 ^m 2 c 2 + (V#) 2 ) = f d d r [m 2 cHp + 9H) . (3.40) 

A second class of invariances involves a reparametrization of the spatial variable, implic- 
itly defined by 

-^ / x -.„, ^\ tan mcuj _ (\ — cos mcui\ .„ i<s 

r -> R(t, r) = r - Q9(t, R) + QCj-r[ (3.41) 

mc \ cos moo J 

„ , . 9(t, R) — mcu) ■ r sinmcu; 

9(t, r -> 9„(t, r = -^-^ 3.42 

cosmcw 



Here a; is a vectorial parameter, O = u)/uj,uj = y/uj 2 . The time-independent generator of 
the infinitesimal transformation reads 

G = J d d r (m 2 c 2 vp + 9pV9) = f d d r (m 2 c 2 rp + 9V) . (3.43) 

With the addition of D and G to the previous generators, the Poincare algebra in (rf+1, 1) 
dimension is reconstructed, and the transformation laws (|3.38| ), (|3.39| ), (|3.41| ), (|3.42| ) ensure 



that (t, r, 9) transforms as a (d+2)-dimensional Lorentz vector (in Cartesian components) [0. 
Note that this symmetry also holds in the free, a = 0, theory. 

Because of the extended symmetry, one can generate new solutions from old ones since 



both 9 W and 9 W in ( |3.39| ) and ( |3.42| ) solve the equation of motion if 9 does. 

A remarkable fact is that the nonrelativistic limit of the above relativistic and interacting 
model precisely corresponds to the nonrelativistic interacting model discussed previously. 
This is easily seen from (|3.30| ), which in the limit gives (|3.1| ), with the help of (|2.43|) . 
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Lorentz 



d 

d£ 



/ d d r mchp + / d d r 



OnrP - P 



'V9 



NR, 



2m 



2mp 



-mc 2 N + I$jgZ 



Equivalently, when p is eliminated, we have from ( |3.36|) and ( [2.43 ) 
U^ -a I d d r \/2m~9 



?NR 



[V9 



NR) 



J X=a 2 /2m 



(3.44) 



(3.45) 



Correspondingly, the equation of motion ( |3.37|) goes over into (|3.17|) . 

It is easy to exhibit solutions of the relativistic theory, which reduce to solutions of the 
nonrelativistic equations that were given previously. The following profiles solve ( |3.37| ). 



G(t,r) = -mc^JcH 2 + r 2 /{d - 1) . 
With ( |2.43| ), this reduces to ( (3.X8|) . In one dimension we have 



9{t, x) = —mcycH 2 + cosh 2 kx/k 2 
reducing to (|3.19aj) . The relativistic analog of the lineal solution (|3.22|) is 



6(t, r) = 9(n • r) + mu • r — mctJc 2 + u 2 — (n • u) 2 . 



(3.46) 



(3.47) 



(3.48) 



Note that the above profiles continue to solve (|3.37|) , even when the sign of the square root 
is reversed; but then they no longer possess a nonrelativistic limit. 

Additionally, there exists an essentially relativistic, chiral solution describing massless 
propagation in one direction: 9 can satisfy the wave equation 



when 



U9 = 



(d^9) 2 = constant 



as, for example, with plane waves 

9{t,r) = /(n-r±ct) 
where (<9 M #) 2 vanishes. Then p reads from (|3.35|) 



P = T- 



mc z 



;f 



(3.49a) 
(3.49b) 

(3.50) 
(3.51) 



IV. RELATION TO NAMBU-GOTO ACTION 

The Nambu-Goto action for a rf-brane in (d + 1) spatial dimensions, moving in time on 
(d + 1,1) Minkowski space is 



In- 



N-G 



f d<fd<p l ■ ■ ■ d<f) d VG 



(4.1) 



where G is (— l) d times the determinant of the induced metric 
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dX» OX. . . 

G *P = l^ldT- ( 4 -2) 



Here Greek letters, from the beginning of the alphabet, label the quantities <p a = (<f)°, 0), with 
which the d-brane coordinate X^ is parametrized; 0° is the evolution parameter, and <fi — 
{4> a , a = 1, . . . , d} are the fixed-time, spatial parameters. These d-brane coordinates carry a 
Greek-letter index from the middle of the alphabet, with value for the temporal coordinate 
X° and m for the <i-brane's d + 1 spatial coordinates. X = {X m , m = 1, . . . , d, d + 1}. 

The action is invariant against reparametrizations of the Q , and we make the 
parametrization choice that the d coordinates X m , m = 1, . . . , d, are given by <ft m , which 
we rename r m . For the remaining parameters we use one of two options, "light-cone" and 
"Cartesian" . 

In the light-cone option, we define 

X ± = -^=(X° ± X d+l ) (4.3) 

and for X + choose the parametrization X + = v2Am</>°; also we rename X + as t. The 
remaining coordinate X~ , a function of (fi° = t/y/2\m and = r, is renamed 9(t,r)/m. 
Upon evaluating the determinant G, we see that the Nambu-Goto action coincides with the 
action for ( |3.16|) [I^|. This identity also explains the higher symmetry noted in Eqs. ( |3.10| ) 



and ( p,13| ). Our choice of parametrization does not interfere with invariance against the 
(d + 1, 1) Poincare group, which acts on X^. In the chosen parametrization, the Poincare 
transformation acts nonlinearly, mixing coordinates (£, r) with the field 9. (However, the 
higher symmetry is also enjoyed by the noninteracting theory, A = 0, which is not equivalent 
to the Nambu-Goto model.) 

For the second, Cartesian option the chosen parametrization permits writing X°, which 
is renamed ct, as amc<f) , while the last coordinate, X d+1 , which is a function of 0° = t/am 
and <f> = r, is called 9(t,r)/mc. With these choices the Nambu-Goto action coincides with 
that for the Born-Infeld Lagrangian, Eq. (|3.36 ). Again the higher dynamical symmetry, 



described by Eqs. ( ^.38| )-( p.43| ), is now understood as the covariance of the Nambu-Goto 
variables X M against (d + 1, l)-dimensional Poincare transformations. (But once again, the 
similar invariance of the free theory, a = 0, cannot be related to properties of a Nambu-Goto 
action.) 

Since both the (d, l)-dimensional Galileo-invariant Chaplygin gas equations and 
the (d, l)-dimensional Poincare- invariant Born-Infeld equations correspond to different 
parametrizations of the (d + 1,1) Nambu-Goto action, there must be a transformation - 
recognized as a reparametrization — that takes solutions of one model into the other. This 
transformation may be formulated as follows. 

Given a solution 0NR{t, r) to (|3.17|), we solve for T(t, r) from the equation 



y/2 V mc 2 

Then a solution 0r(£, r) to ( |3.37|) , is given by 



^=(T(t,r) + -L^(T(t,r),r)) =t. (4.4) 



-4(T(t,r) - -^9 NR (T(t,r),r)) = -^9 R {t,r) (4.5a) 

a/9 \ m.r z V / / m.r z 



a/2 V mc z v ' ) mc* 

or 
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6 R (t,r)=mc 2 (V2T(t,r)-t) . 



(4.5b) 



Indeed this mapping produces solutions ( |3.46|) , ( |3.47|) and (|3.48|) from (|3.18|) , (|3.19a|) and 
( |3.22j ) respectively. [In fact both signs of the square root are obtained; also ( |3.48| ) emerges 



from ( 3.22 ) only after a redefinition of and u.] 

Oppositely, given as solution 0r(£, r) to the relativistic Born-Infeld equation (|3.37| ), we 
reparametrize by solving for T(t,r) from 



and then find 0NR(t, r) from 



-L(T( t ,r)-^ fi (T(t, 



r ,r 



mr 



;9 NR (t, r) 



or 



9 NR (t,T) = mc 2 (V2T(t,r) - t 



(4.6) 



(4.7a) 



(4.7b) 



The two transformations are collected in the statement 

e NR (t,r)- -' 



mc 2 (V2T - t) 



9r(T,t) 



mc 



2t-T) 



(4- 



with the instruction that obtaining an expression for 9n R in terms of Or, or vice-versa, 
requires that one of T or t be eliminated in favor of the other. 

The interrelationships may be summarized by the following diagram: 



^Para 

N; 
^^in 



Parameterization-invariant 
Nambu-Goto d-brane 
(d + 1, 1) spacetime 



7 ^ 



iriant\ 
ine 



Light-cone 
parameterization 



/ 



Cartesian 
parameterization 



' Galileo-invariant 



Chaplygin gas 
in (d, 1) spacetime 



Nonrelativistic limit 
< 



V 



Exact transformation 



Poincare-invariant 
Born-Infeld model 
in (d, 1) spacetime 



e 



It is striking that there exists a two-fold relationship between the Chaplygin gas and the 
Born-Infeld model. First, there is the exact mapping, given in ( f4.8| ), of one into the other. 
Second, the nonrelativistic limit (|3.44|) , (|3.45|) of the latter produces the former. 
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V. ONE-DIMENSIONAL MOTION 

In one spatial dimension, the motion of these systems simplifies and they become integrable. 
We shall give a self-contained demonstration of integrability and derive the integrals of mo- 
tion in a compact form, stressing the connection between the relativistic and nonrelativistic 
case. 

In one dimension the requirement that the local momentum is irrotational poses no 
restriction. We shall use as phase space variables the local momentum p(t, x) and the particle 
density p(t,x). The equal-time Poisson bracket fl2.25|) becomes 



{p(x), p(x')} = S'(x — x') . (5.1) 

Note that p and p in this relation are on an equal footing and are governed by a nonlocal 
canonical 1-form | J dxdyp(x)e(x — y)p(y)dt, where e is the antisymmetric step function. 
We shall consider local integrals, that is, quantities of the form 



dxT{p(x), p(x)) (5.2) 

with T a local function of p and p. The Poisson bracket of two such integrals F and G is 
calculated through (|5.1|) as 

{F, G} = - I dx {TppQp + TppQ p )d x p + [T PV Q P + T w Q p )d x p (5.3) 

where we suppressed the dependence on x and subscripts indicate partial derivative. If 
the above integrand is a total x-derivative then (with appropriate boundary conditions) the 
integral will vanish and F and G will be in involution. For this we need the curl-free condition 

ly ppyp T J" ppV p)p \y ppi/p T J~ ppSj p) p V / 

or, finally 

r a 

(5.5) 



J~ pp __ ypp 



j~pp ypp 
Choosing one of the integrals to be the Hamiltonian H = J dxTi, the well-known condition 

^p = Hjv (5 6) 

J~pp iLpp 

guarantees that F is a constant of the motion. If we recover a set of such integrals satisfying 
( |5.6| ) then they will obviously also satisfy ( |5.5|) among themselves. Therefore, constants of 
motion will automatically be in involution. 

For the nonrelativistic case the Hamiltonian density ( |3.2j ), ( p.4[ ) is 

n = pf- + - (5.7) 

2m p 
and therefore the integrals of motion are generated by functions that satisfy 

p^Tpp = 2\mT vv . (5.8) 



This can readily be solved by separation of variables. We prefer, however, simply to give its 
general solution in terms of two arbitrary functions / and g of one variable: 



_ „/ \ / 2\m\ / \/2Am\ , N 

F = pf[p+ X )+pg[p- x ) • (5-9) 

We essentially get two infinite towers of integrals. Choosing, e.g., f(z) 
g(z) = z n , f = we get the integrals 



i±=f*xp(p±^y. (s.io) 



All the integrals presented in ||, |16| can be identified as linear combinations of the 1^- As 
stated, the I* are all in involution and demonstrate the complete integrability of the system. 
The Hamiltonian, in particular, is included as 4H = l£ + 1% and the total momentum as 
2P = J 1 + + /f. 
The quantities 



R ± =p±^™ (5.11) 

P 

appearing above are known as Riemann coordinates. The equations of motion for this system 
(continuity and Euler) are summarized in terms of R±: 

R± = R= R'+ . (5.12) 

m 

Although this formulation for the equations is known, the relation to the constants of motion 
does not seem to appear in the literature. 

Another fluid system for which the equations of motion, expressed in terms of Riemann 
coordinates, take a form similar to ( |5.12j ) possesses a potential cubic in p: V(p) = £p 3 /3. 
This also arises in a collective, semiclassical description of nonrelativistic free fermions, where 
the cubic potential reproduces fermion repulsion [|17| |. In this case, the Riemann corrdinates 
read 



R±=p±V2£mp (5.13) 

and, in contrast to ( |5.12| ), they decouple in the the equations of motion: 

R± = R±R' ± . (5.14) 

m 

Indeed it is seen that R± satisfy essentially the free Euler equation [compare ( |1.2| ) and 
identify R± with v]. Consequently (|5.14|) is solved by analogs of Q2.2|) - ([2.4|) . 

Both these examples are special limiting cases of a more general system, with potential 

V(p) = ^P±M + Ap + B . (5.15) 

{p + af 

(The terms A, B correspond to a dynamically trivial part that does not alter the equations 
of motion.) The Chaplygin gas corresponds to a = A = B = 0, while the cubic potential is 
recovered in the limit 
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A = i(£ A = Ua 2 B = -Ua 3 , a -> oo . (5.16) 



3*"* ~ 3 

The Riemann coordinates are 



i?±=p±^— - 5.17 

p + a 



and the conserved integrals of this system are given by functions of R±: 



*• t w/ V2Am\ . , / V2Am\ /K< . 

V p + a J \ p + a ) 



The physical meaning of this general system is not clear. 

We conclude the discussion on the one- dimensional nonrelativistic Chaplygin gas by 
presenting a set of moving solutions, which are the Galileo boosts of the static solutions 
( gg| ), ( ET23D in Section III. These read 

\/2\m 

p = p(x — ut) , p = -. r (5.19) 

\p — mu\ 

with p(x — ut) an arbitrary function of x — ut (provided it never equals mu). Clearly this is 
a constant-profile solution moving with a velocity u. 

For the relativistic Born-Infeld system, the Hamiltonian density (|3.29|) is 



H = ^p^c 2 + a 2 ^fm 2 c 2 + p 2 . (5.20) 

Relation ( |5.6|) for the conserved integrals reads 

m 2 (p 2 c 2 + a 2 ) 2 T pp = a 2 (m 2 c 2 + p 2 ) 2 T vv . (5.21) 

This can be solved by separation of variables. We prefer again, however, to define Riemann 
coordinates and give the solution in terms of arbitrary functions of one variable, just as in 
the nonrelativistic case. The relevant combinations here are 

pc p 

R± = 4>n ±4> P < \>p = arctan — <fi p = arctan — (5.22) 

a mc 

and the equations of motion read, in terms of these, 

R± = T c(cosR T )R' ± (5.23) 

while the general solution to ( |5.21| ) is 



T 



mca 



cos <p p cos <p p 



/(</V + p )+2(0p-0p) • (5-24) 



By choosing, as before, simple monomials for / and g we get the tower of conserved quantities 

It = mca [ dx (0P ± 0p) " . (5.25) 

J COS (fi p COS (ftp 

The Hamiltonian is included as H = I Q . The momentum, on the other hand, is an infinite 
series in the above integrals, requiring / and g to be exponential functions. We can give 
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an alternative tower of complex conserved integrals involving only algebraic functions of p 
and p: 



I* = / dx \J p 2 c 2 + a 2 ^Jm 2 c 2 +p 2 (pc ± ia) n (p + imc) n . (5.26) 

Then the Hamiltonian is H = I while the total momentum, total particle number and 
integral of the local momentum are contained in the real and imaginary parts of If and I{~ 
(as well as a trivial constant). It can be checked that the integrals (|5.26j ) go over to the 
nonrelativistic ones ( |5.10| ) in the limit c — > oo upon proper rescaling. 

In the relativistic model p need not be constrained to be positive (negative p could be 
interpreted as antiparticle density). The transformation p — > —p, p — > —p is a symmetry 
and can be interpreted as charge conjuguation. Further, p and p appear in an equivalent 
way. As a result, this theory enjoys a duality transformation: 

a , mc2 tr nw\ 

p -► ±—p P -► ± P 5.27 

mc z a 

Under the above, both the canonical structure and the Hamiltonian remain invariant. Solu- 
tions are mapped in general to new solutions. Note that the nonrelativistic limit is mapped 
to the ultra- relativistic one under duality. Self-dual solutions p = ±-%p satisfy 

p = Tcp' (5.28) 

and are, therefore, the chiral relativistic solutions that were presented at the end of Sec- 
tion III. In the self-dual case, when p is eliminated from the canonical 1-form and from the 
Hamiltonian, one arrives at an action for p, which coincides (apart from irrelevant constants) 



with the self-dual action, constructed some time ago |1| : 



p=^p 



2mc 2 



lJd X dyp( X )e( X -y)p(y)dt-Jd X ^J^ 2 ^^7 2 dt] 

\Jdxdyp(x)e{x-y)p(y)dt-^j 'dx (p 2 (x) + ^) dtj . (5.29) 

A set of constant-profile solutions can be found by Lorentz-boosting the static solution 
( |3.48 ). Their most general form is 



a up± mcVc 2 - u 2 
p = p(x-ut), p= (5.30) 

c pv c — u T rncu 

with p( X — ut), again, a general function. Note that the two choices of sign in the above 
formula are related by charge conjuguation. In the extreme relativistic case u — > c this 
solution goes over to the chiral relativistic solution ( [3.51 ). The set of these solutions is 
closed under duality transformations. 

We shall conclude by presenting an explicit mapping between the relativistic and the 
nonrelativistic theories in one dimension, which demonstrates their kinematical equivalence. 
Note that the relativistic equations of motion ( 5.23| ) become, in terms of cosi?±, 



d t cosR± = ^fc(cosR :p )d x cosR± . (5.31) 
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These are essentially identical to the equations of motion for the nonrelativistic Riemann 
coordinates ( |5.12| ). In fact, putting 

R± = ±mc cos R± (5.32) 

we see that the R± obey the nonrelativistic equations ( |5.12| ). Expressing R± and R± in terms 
of the corresponding nonrelativistic and relativistic variables produces a mapping between 
the two sets. Call Pnr, Pnr and pr, pr the local momentum and density of the nonrelativistic 
and the relativistic theory, respectively. Then the mapping is 

V2Xm 2 PRPR /k qqx 

Pnr = ^rl R , p NR = mc —— (5.33) 

am z c z Hr 



where TCr = Jp 2 R c 2 + a 2 Jm 2 c 2 + p R is the relativistic Hamiltonian density. As can be 
checked by direct algebra, this maps the relativistic equations of motion to the nonrelativistic 
ones. Since the combinations of p R and pR that appear in ( |5.33|) are duality- invariant, the 
mapping of solutions is two-to-one. Note that the constant-profile relativistic solutions ( |5.30|) 
are mapped to the corresponding nonrelativistic ones ( |5.19 



We stress that the transformation ( [5.33J ) is not canonical, since it does not preserve the 



Poisson brackets. Accordingly, it does not map the relativistic Hamiltonian into the nonrel- 
ativistic one. This is a manifestation of the bi-Hamiltonian structure of these systems, since 
there are now two pairs of Hamiltonian and canonical structure (the standard nonrelativis- 
tic one and the one obtained through this mapping) that lead to the same nonrelativistic 
equations of motion. 

We note that a similar mapping between the nonrelativistic system and the relativistic 
one in light-cone coordinates was presented previously by Verosky [[TJJ. The mapping (]5T! 



then, can be considered as the analog of Verosky's transformation for the Lorentz system, 
although it cannot be obtained from it in any straightforward way. 

The existence of an infinite set of constants of motion for both the nonrelativistic Chap- 
lygin gas, ( ^9|) and (|5.10|) , as well as for the relativistic Born-Infeld model, fl5.24|) -( [5.26l ), 
signals the complete integrability of these theories. The actual integration of the equations of 
motion can be carried out only indirectly. First of all, since both these d = 1 models descend, 
via alternate parametrization choices, from the Nambu-Goto string (1-brane) on the plane 
(2-space), the explicit integration of the latter j|] allows presenting solutions of the former 
two in terms of two arbitrary functions. Alternatively, the Chaplygin gas equations can be 
combined, after a Legendre transformation, into a linear, second-order partial differential 
equation ||, whose general solution, in terms of two arbitrary functions, is known explic- 



itly p0[ . The Born-Infeld solution can then be constructed with the help of transformations 



described in Sections IV and V. 



22 



REFERENCES 

[1] J. Goldstone, (unpublished). 

[2] M. Bordemann and J. Hoppe, Phys. Lett. B317, 315 (1993); B325, 359 (1994). 

[3] Y. Nutku, J. Math. Phys. 28, 2579 (1987); P. Olver and Y. Nutku, J. Math. Phys. 29, 

1610 (1988); M. Arik, F. Neyzi, Y. Nutku, P. Olver and J. Verosky, J. Math. Phys. 30, 

1338 (1989). 
[4] See e.g. D. Bazeia, Phys. Rev. D (in print). 
[5] J. Brunelli and A. Das, Phys. Lett. A235, 597 (1997). 
[6] C. Eckart, Phys. Rev. 54, 920 (1938); recent work: A. Schakel, Mod. Phys. Lett. BIO, 

999 (1996); N. Ogawa, preprint, |hep-th/980lTT5 . 



[7] L.D. Faddeev and R. Jackiw, Phys. Rev. Lett. 60, 1692 (1988). 

[8] In conventional treatments of relativistic hydrodynamics, as in S. Weinberg, Gravitation 

and Cosmology (Wiley, New York, 1972), p is a scalar rather than the time- component 

of a Lorentz vector. This distinction makes no difference in the noninteracting case, but 

our choice allows introducing self-interactions in a Lorentz-invariant manner. 
[9] L. Landau and E. Lifschitz, Fluid Mechanics , 2nd ed. (Pergamon, Oxford UK, 1987). 
[10] E. Madelung, Z. Phys. 40, 322 (1926); recent work: E. Merzbacher, Quantum Mechanics, 

3rd ed. (Wiley, New York, 1998). 
[11] D. Bazeia and R. Jackiw, Ann. Phys. (NY), 270, 246 (1998); R. Jackiw and A. Poly- 

chronakos, Faddeev Festschrift, preprint, |hep-th/9809123| . 
[12] A. Jevicki, Phys. Rev. D 57, 5955 (1998). 
[13] Bordemann and Hoppe, second paper in Ref. [2]. 
[14] This result was established with a hodographic transformation for d = 2 by Bordemann 

and Hoppe (first paper in Ref. [2]) and for arbitrary d by Jackiw and Polychronakos, 

Ref. [11]; see also Bazeia, Ref. [4]. The present, simple approach was shown to us by 

Hoppe; see also J. Hoppe, Phys. Lett. B329, 10 (1994). 
[15] This result was established with a hodographic transformation for d = 2 by Bordemann 

and Hoppe (second paper in Ref. [2]); see also Hoppe, Ref. [14]. 
[16] Nutku, Olver and Nutku, Ref. [3]. 
[17] A. Jevicki and B. Sakita, Phys. Rev. D22, 467 (1980) and Nucl. Phys. B165, 511 (1980); 

J. Polchinski, Nucl. Phys. B362, 25 (1991); J. Avan and A. Jevicki, Phys. Lett. B266, 

35 (1991) and B272, 17 (1991). 
[18] R. Floreanini and R. Jackiw, Phys. Rev. Lett. 59, 1873 (1987). 
[19] As cited in Olver and Nutku, Ref. @. 



[20] Jackiw and Polychronakos, Ref. |TT 



23 



